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Abstract 

Silver  nanoparticles  (AgNPs)  are  widely  produced  and  used.  Because  of  their 
potential  toxicity  and  the  possibility  of  their  release  into  the  environment,  it  is  important 
to  understand  the  factors  involved  in  their  transport;  particularly  how  they  may  move  in 
groundwater,  which  is  a  potential  pathway  to  human  and  environmental  receptors.  By 
passing  a  solution  of  15  mg/L  AgNPs  with  an  average  size  of  17  mn  through  sand-filled 
glass  columns,  this  study  looks  at  the  physical  effects  of  flow  rate  and  media  size  on 
transport.  Three  different  sand  sizes  (<  0.074,  0.21-0.297,  and  0.4-2. Onnn)  were  used,  as 
well  as  two  different  flow  rates  (1  and  4  mL  min-1).  Results  showed  that  flow  rate  had 
little  effect  on  the  transport  of  AgNPs.  As  had  been  found  in  previous  studies,  transport 
through  fine  sand  (<  0.074  mm)  was  inhibited,  and  no  AgNPs  were  observed  in  the 
column  effluent.  In  contrast  to  other  studies  and  conventional  filtration  theory,  however, 
it  was  observed  that  more  AgNPs  were  captured  in  the  coarse  sand  (0. 4-2.0  mm)  than  the 
medium  sand  (0.21-0.297  mm).  This  deviation  from  conventional  filtration  theory  was 
attributed  to  the  more  heterogeneous  size  distribution  of  the  coarse  sand  compared  to  the 
medium  sand. 
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INFLUENCE  OF  MEDIA  SIZE  AND  FLOW  RATE  ON  THE  TRANSPORT  OF 
SILVER  NANOPARTICLES  IN  SATURATED  POROUS  MEDIA:  LABORATORY 

EXPERIMENTS  AND  MODELING 

I.  Introduction 

Research  on  silver  nanoparticles  (AgNPs)  has  increased  recently  because  of  their 
extensive  use  in  consumer  products  (Benn  et  al.  2010;  Impellitteri  et  al.  2009),  medical 
and  pharmaceutical  supplies  (Faunce  and  Watal  2010),  environmental  applications  and 
electronics  (Kim  et  al.  2007;  Tolaymat  et  al.  2010;  Wijnhoven  et  al.  2009).  Some  of 
these  studies  have  shown  that  AgNPs  may  have  toxic  effects  on  the  environment,  aquatic 
organisms  and  human  health  (Kim  et  al.  2008;  Wijnhoven  et  al.  2009;  Schrand  et  al. 
2010;  Kennedy  et  al.  2010;  Shaw  and  Handy  2011).  Because  of  their  toxicity  and 
ubiquity,  it  is  necessary  to  better  understand  the  fate  and  transport  of  AgNPs  in  the 
environment.  In  particular,  understanding  the  fate  and  transport  of  AgNPs  in 
groundwater  is  necessary,  as  groundwater  is  an  important  pathway  of  AgNPs  to  human 
and  environmental  receptors  (Flory  et  al.,  2013). 

Filtration  theory  presented  by  Yao  et  al.  (1971)  posits  that  there  are  three  main 
mechanisms  that  affect  the  capture  of  particles  in  porous  media.  The  three  mechanisms 
are  advection  (interception),  diffusion  and  gravity  (sedimentation).  Yao  et  al.  (1971) 
state  that  removal  of  particles  in  the  size  range  of  AgNPs  (<100  nm)  would  be  controlled 
by  diffusion. 

Because  diffusion  is  the  primary  mechanism  controlling  AgNP  transport,  as  the 
velocity  of  the  fluid  increases  the  percentage  of  particles  or  mass  that  the  filter  captures 
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decreases.  In  addition,  as  the  media  size  decreases  the  efficiency  of  the  filter  should 
increase.  This  means  that  for  AgNP  transport  through  porous  media,  more  particles  will 
be  collected  as  the  media  diameter  gets  smaller. 

Li  et  al.  (2008)  conducted  studies  on  the  effect  that  two  flow  rates  had  on  the 
transport  of  Fullerene  nanoparticles  with  four  fractions  of  Ottawa  sand.  In  accordance 
with  filtration  theory,  they  found  that  as  the  media  size  or  flow  rate  decreased,  a  higher 
fraction  of  mass  was  retained  in  the  media. 

Multi-walled  carbon  nanotube  (MWCNT)  transport  was  studied  by  Mattison  et  al. 
(2011).  They  found  with  different  flow  rates  and  sand  sizes  that  MWCNT  retardation 
increased  with  decreasing  collector  size,  which  is  in  accord  with  colloid  filtration  theory. 
In  their  study  two  pulses  of  MWCNTs  were  injected  into  the  column.  They  found  that 
deposition  rate  decreases  as  their  experiments  progressed  and  this  could  not  be  predicted 
with  traditional  colloid  filtration  theory.  They  indicated  that  chemical  and  physical 
heterogeneity,  both  of  which  would  be  expected  in  subsurface  environments,  could  affect 
transport  and  are  important  topics  for  future  research. 

Some  studies  have  been  conducted  to  detennine  the  effect  of  media  heterogeneity 
on  the  transport  of  various  bacteria,  viruses  and  colloids.  Silliman  et  al.  (2001)  found  that 
some  strains  of  bacteria  traveled  through  a  homogenous  porous  media  better  than  a  more 
heterogeneous  one.  One  important  finding  was  that  as  media  heterogeneity  increased,  the 
percentage  of  mass  leaving  their  column  decreased.  In  their  study  they  quantified  their 
heterogeneity  by  the  hydraulic  conductivity  and  range  of  grain  sizes.  Also  Redman  et  al. 
(2001)  evaluated  the  impact  of  media  heterogeneity  on  pathogen  filtration.  They  found 
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that  heterogeneity  at  the  scale  of  the  pathogen  or  the  collector  (microscale)  leads  to  a 
slow  decay  of  concentration  versus  the  fast  exponential  decay  predicted  by  standard 
filtration  theory.  They  also  found  that  heterogeneity  at  the  filter  scale  (macroscale)  can 
have  large  effects  on  the  transport  of  viruses.  This  macroscale  heterogeneity  describes 
the  difference  in  size  of  grains  in  the  filter  versus  the  microscale  which  describes  how 
smooth  or  rough  individual  particles  are  within  the  media. 

Previous  work  that  has  investigated  the  effect  of  flow  rate  and  media  size  on 
engineered  nanopartices  (ENPs)  transport  has  used  various  ENP  types.  This  study 
focuses  on  how  flow  rate  and  media  grain  size  affects  the  transport  of  AgNPs,  which  due 
to  their  toxicity  and  ubiquity  are  of  particular  importance.  Specifically,  this  study  will 
focus  on  what  physical  processes  are  relevant  to  the  transport  of  AgNPs  in  porous  media 
and  what  model  governing  equations  and  parameter  values  are  appropriate  to  describe 
AgNP  fate  and  transport  in  porous  media  under  varying  hydrogeological  (flow  rate  and 
media  size)  conditions. 

Problem  Statement 

While  several  studies  have  looked  at  ENP  transport  in  porous  media  at  different 
media  grain  sizes  and  flow  rates  none,  to  the  authors’  knowledge,  have  conducted  a 
focused  investigation  on  how  these  two  factors  in  particular  affect  transport  of  AgNPs. 
This  study  focuses  on  how  flow  rate  and  media  grain  size  affects  the  transport  of  AgNPs. 
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Research  Objectives 

The  objective  of  this  study  was  to  further  understand  how  AgNPs  are  transported 
in  groundwater.  Specifically,  the  following  questions  are  answered: 

1 .  What  physical  processes  are  relevant  to  transport  of  AgNPs  in  sand? 

2.  What  model  governing  equations  and  parameter  values  are  appropriate  to  describe 
AgNP  fate  and  transport  in  porous  media  under  varying  hydrogeological  (flow  rate  and 
media  size)  conditions? 

Preview 

This  thesis  is  written  using  the  scholarly  article  format.  Chapter  II  contains  a 
manuscript  intended  for  submission  to  a  peer-reviewed  journal.  The  manuscript  includes 
an  abstract,  along  with  introduction,  materials  and  methods,  and  results  and  discussion 
sections.  Chapter  III  provides  a  summary  of  the  primary  findings  discussed  in  the  article, 
as  well  as  related  topics  for  future  study. 
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II.  Scholarly  Article 


Influence  of  media  size  and  flow  rate  on  the  transport  of  silver  nanoparticles  in 
saturated  porous  media:  Laboratory  experiments  and  modeling 

Travis  Meidinger  Sushil  R.  Kanel,  LeeAnn  Racz,  Chelsea  Marcum,  Mark  N.  Goltz 

Abstract 

Silver  nanoparticles  (AgNPs)  are  widely  produced  and  used.  Because  of  their 
potential  toxicity  and  the  possibility  of  their  release  into  the  environment,  it  is  important 
to  understand  the  factors  involved  in  their  transport;  particularly  how  they  may  move  in 
groundwater,  which  is  a  potential  pathway  to  human  and  environmental  receptors.  By 
passing  a  solution  of  15  mg/L  AgNPs  with  an  average  size  of  17  mn  through  sand-filled 
glass  columns,  this  study  looks  at  the  physical  effects  of  flow  rate  and  media  size  on 
transport.  Three  different  sand  sizes  (<  0.074,  0.21-0.297,  and  0.4-2. 0mm)  were  used,  as 
well  as  two  different  flow  rates  (1  and  4  mL  min-1).  Results  showed  that  flow  rate  had 
little  effect  on  the  transport  of  AgNPs.  As  had  been  found  in  previous  studies,  transport 
through  fine  sand  (<  0.074  mm)  was  inhibited,  and  no  AgNPs  were  observed  in  the 
column  effluent.  In  contrast  to  other  studies  and  conventional  filtration  theory,  however, 
it  was  observed  that  more  AgNPs  were  captured  in  the  coarse  sand  (0. 4-2.0  mm)  than  the 
medium  sand  (0.21-0.297  mm).  This  deviation  from  conventional  filtration  theory  was 
attributed  to  the  more  heterogeneous  size  distribution  of  the  coarse  sand  compared  to  the 
medium  sand. 
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Introduction 


Silver  nanoparticles  (AgNPs)  have  been  studied  much  as  of  late  because  of  their 
extensive  use  in  consumer  products  (Benn  et  al.  2010;  Impellitteri  et  al.  2009),  medical 
and  phannaceutical  supplies  (Faunce  and  Watal  2010),  environmental  applications  and 
electronics  (Kim  et  al.  2007;  Tolaymat  et  al.  2010;  Wijnhoven  et  al.  2009).  Some  of 
these  studies  have  shown  that  AgNPs  have  potentially  toxic  effects  for  the  environment, 
fish  and  humans  (Kim  et  al.  2008;  Wijnhoven  et  al.  2009;  Schrand  et  al.  2010;  Kennedy 
et  al.  2010;  Shaw  and  Handy  2011).  Because  of  their  toxicity  and  ubiquity,  it  is 
necessary  to  better  understand  the  fate  and  transport  of  AgNPs  in  the  environment.  In 
particular,  understanding  the  fate  and  transport  of  AgNPs  in  groundwater  is  necessary,  as 
groundwater  is  an  important  pathway  of  AgNPs  to  human  and  environmental  receptors 
(Flory  et  al.,  2013). 

Filtration  theory  presented  by  Yao  et  al.  (1971)  posits  that  there  are  three  main 
mechanisms  that  affect  the  capture  of  particles  in  porous  media.  The  three  mechanisms 
are  advection  (interception),  diffusion  and  gravity  (sedimentation).  Yao  et  al.  (1971) 
state  that  removal  of  particles  in  the  size  range  of  AgNPs  (<100  nm)  would  be  controlled 
by  diffusion. 

Because  diffusion  is  the  primary  mechanism  controlling  AgNP  transport,  as  the 
velocity  of  the  fluid  increases  the  percentage  of  particles  or  mass  that  the  filter  captures 
decreases.  In  addition,  as  the  media  size  decreases  the  efficiency  of  the  filter  should 
increase.  This  means  that  for  AgNP  transport  through  porous  media,  more  particles  will 
be  collected  as  the  media  diameter  gets  smaller. 
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Li  et  al.  (2008)  conducted  studies  on  the  effect  that  two  flow  rates  had  on  the 
transport  of  Fullerene  nanoparticles  with  four  fractions  of  Ottawa  sand.  In  accordance 
with  filtration  theory,  they  found  that  as  the  media  size  or  flow  rate  decreased,  a  higher 
fraction  of  mass  was  retained  in  the  media. 

Multi-walled  carbon  nanotube  (MWCNT)  transport  was  studied  by  Mattison  et  al. 
(2011).  They  found  with  different  flow  rates  and  sand  sizes  that  MWCNT  retardation 
increased  with  decreasing  collector  size,  which  is  in  accord  with  colloid  filtration  theory. 
In  their  study  two  pulses  of  MWCNTs  were  injected  into  the  column.  They  found  that 
deposition  rate  decreases  as  their  experiments  progressed  and  this  could  not  be  predicted 
with  traditional  colloid  filtration  theory.  They  indicated  that  chemical  and  physical 
heterogeneity,  both  of  which  would  be  expected  in  subsurface  environments,  could  affect 
transport  and  are  important  topics  for  future  research. 

Some  studies  have  been  conducted  to  detennine  the  effect  of  media  heterogeneity 
on  the  transport  of  various  bacteria,  viruses  and  colloids.  Silliman  et  al.  (2001)  found  that 
some  strains  of  bacteria  traveled  through  a  homogenous  porous  media  better  than  a  more 
heterogeneous  one.  One  important  finding  was  that  as  media  heterogeneity  increased,  the 
percentage  of  mass  leaving  their  column  decreased.  In  their  study  they  quantified  their 
heterogeneity  by  the  hydraulic  conductivity  and  range  of  grain  sizes.  Also  Redman  et  al. 
(2001)  evaluated  the  impact  of  media  heterogeneity  on  pathogen  filtration.  They  found 
that  heterogeneity  at  the  scale  of  the  pathogen  or  the  collector  (microscale)  leads  to  a 
slow  decay  of  concentration  versus  exponential  decay  predicted  by  standard  filtration 
theory.  They  also  found  that  heterogeneity  at  the  filter  scale  (macroscale)  can  have  large 
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effects  on  the  transport  of  viruses.  This  macroscale  heterogeneity  describes  the 
difference  in  size  of  grains  in  the  filter  versus  the  microscale  which  describes  how 
smooth  or  rough  an  individual  particles  are  within  the  media. 

Previous  work  that  has  investigated  the  effect  of  flow  rate  and  media  size  on 
engineered  nanoparticles  (ENPs)  transport  has  used  various  ENP  types.  This  study 
focuses  on  how  flow  rate  and  media  grain  size  affects  the  transport  of  AgNPs,  which  due 
to  their  toxicity  and  ubiquity  are  of  particular  importance.  Specifically,  this  study  will 
focus  on  what  physical  processes  are  relevant  to  the  transport  of  AgNPs  in  porous  media 
and  what  model  governing  equations  and  parameter  values  are  appropriate  to  describe 
AgNPs  fate  and  transport  in  porous  media  under  varying  hydrogeological  (flow  rate  and 
media  size)  conditions. 

Materials  and  Methods 

Materials 

Two  types  of  sand  were  purchased  for  these  experiments.  The  first  was  50-70 
mesh  sand  (0.297-0.21  mm)  from  Sigma- Aldrich  (St.  Louis,  MO).  The  second  sand  was 
Quikrete  50  lb.  Commercial  Grade  Medium  Sand  (Atlanta,  GA)  with  an  advertised  size 
range  of  20-50  mesh  (0.841-0.297  mm).  A  1000  mg  L'1  suspension  of  PVP  doped 
AgNPs  (hereafter  referred  to  as  AgNPs)  (99.99%  pure)  was  purchased  from  US  Research 
Nanomaterials,  Inc.  (Houston,  TX).  ACS  reagent  grade  HN03  was  purchased  from 
Sigma-Aldrich  (St.  Louis,  MO).  The  reagents  were  used  as  received. 
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Sand  Preparation 

The  Quikrete  sand  was  sieved  to  produce  a  10-40  mesh  fraction  (0.4-2  mm)  and  a 
>200  mesh  (<0.074  mm)  fraction.  The  two  size  fractions  were  stored  in  individual  sealed 
glass  containers.  Each  size  fraction  was  then  acid  washed.  The  sand  was  put  in  a 
Buechner  funnel  with  Whatman  No  42  filter  paper  on  top  of  the  sand.  The  sand  was  then 
rinsed  with  water  and  the  water  was  vacuumed  off.  The  funnel  was  then  filled  with 
15.9N  HNO3  and  allowed  to  sit  for  5  min.  The  sand  was  then  rinsed  with  aliquots  of 
reverse  osmosis  (RO)  purified  water  until  it  reached  a  pH  of  7.0.  It  was  then  dried  in  an 
oven  at  300°  C  for  two  hours.  Once  dry,  the  sands  were  placed  in  clean  sealed  glass 
containers.  The  Sigma-Aldrich  50-70  mesh  sand  was  used  as  it  came  from  the  supplier, 
without  additional  preparation. 

Column  Experiments 

A  glass  column  (2.5  cm  inner  diameter  x  5  cm  length)  was  packed  with  42g,  45 g 
or  44g  of  10-40  mesh  (0.4-2  mm),  50-70  mesh  (0.21-0.297  mm),  or  >200  mesh  (<0.074 
mm)  sand  respectively  (referred  to  as  coarse,  medium  and  fine  sand  throughout  this 
paper),  and  0.1  g  of  glass  wool  was  added  to  each  end  of  the  column  to  keep  the  media  in 
the  column.  The  column  was  packed  wet  for  each  experiment  to  ensure  no  cross 
contamination  between  experiments.  Deionized  (DI)  water  was  passed  through  the 
column  for  approximately  30  minutes  to  ensure  there  were  no  leaks  or  air  bubbles  in  the 
column  or  tubing.  The  flow  rate  was  set  on  a  peristaltic  pump  (0-100  rpm,  MasterFlex) 
to  1  or  4  mL  min'1  depending  on  the  experiment  being  run.  The  flow  rate  was  measured 
by  measuring  the  volume  of  solution  that  was  captured  by  the  fraction  collector  as  a 
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function  of  time  for  each  experiment.  Once  the  desired  flow  rate  was  reached,  tracer 
solution  (60  mL  of  3  mM  aqueous  KC1  solution)  and/or  AgNP  (60  mL  of  15  mgL'1  AgNP 
suspension)  was  passed  through  the  column.  The  solutions  were  added  for  1 5  minutes  at 
4  mL  min-1  or  60  min  at  1  mL  min-1 .  After  the  allotted  time,  DI  water  was  then  passed 
through  the  column  for  an  additional  30  or  120  minutes  for  the  4  and  1  mL  min'1  runs, 
respectively.  The  flow  rate  was  monitored  throughout  the  experiment  and  adjusted  if 
needed.  The  effluent  was  passed  through  an  open  cell  in  an  Ultraviolet-visible  (UV-vis) 
spectrophotometer  (UV-Cary  60,  Agilent).  Readings  were  taken  on  the  UV-vis  every  30 
seconds  for  the  4  mL  min'1  experiment  and  every  2  minutes  for  the  1  mL  min'1 
experiment.  The  UV-vis  was  zeroed  prior  to  each  run  with  DI  water.  In  addition  a 
calibration  curve  was  developed  by  using  AgNPs  at  a  concentration  of  50,  25,  and  12.5 
mgL'1.  The  limit  of  detection  for  the  UV-vis  is  approximately  0.5  mgL'1  and  a  limitation 
of  the  instrument  is  it  only  detects  silver  particles  and  does  not  measure  total  silver.  All 
column  experiments  were  performed  in  duplicate  and  average  concentration  values  are 
reported. 

A  plot  showing  normalized  effluent  concentration  (the  concentration  at  the 
column  outlet  divided  by  the  influent  concentration)  versus  dimensionless  time  (pore 
volume,  PV)  was  developed  for  each  experiment.  This  dimensionless  breakthrough 
curve  (BTC)  of  effluent  concentration  vs  time  allows  for  easy  comparison  between 
experiments. 

Zeroth  and  first  moments  of  the  dimensionless  breakthrough  curves  were 
calculated.  The  zeroth  moment,  which  is  the  area  under  the  breakthrough  curve, 
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represents  the  total  mass  of  tracer  or  AgNPs  that  passed  through  the  column  and  was  in 
the  effluent.  Comparing  the  zeroth  moment  of  the  AgNP  BTC  with  the  zeroth  moment  of 
the  nonreactive  tracer  allows  us  to  quantify  the  fraction  of  AgNPs  that  was  retained  by 
the  media.  Similarly,  the  first  moment  represents  the  mean  time  tracer  and  AgNPs  were 
in  the  column.  Comparing  the  first  moment  of  the  AgNP  BTC  with  the  first  moment  of 
the  nonreactive  tracer  allows  us  to  quantify  the  retardation  of  the  AgNPs  as  they  are 
transported  through  the  column. 

Analysis 

The  AgNPs  used  in  these  experiments  were  characterized  and  reported  in  Flory  et 
al.  (2013).  High  resolution  transmission  electron  microscopy  (HRTEM)  was  used  to 
measure  the  shape  and  size  of  AgNPs  on  a  Hitachi  H-7600  transmission  electron 
microscope  at  an  accelerating  voltage  of  100  kV.  In  addition,  TEM  images  where  taken 
on  two  samples  that  had  passed  through  the  column.  X-ray  diffraction  (XRD)  (Bruker 
D8  Advanced)  with  a  CuK-alpha  source  was  used  to  identify  any  metal  oxides  present  on 
the  three  sizes  of  sand. 

The  contents  of  the  column  were  retained  and  analyzed  to  measure  the  amount  of 
silver  associated  with  either  the  water  or  the  amount  captured  by  the  sand  that  was  left  in 
the  column.  After  each  experiment  the  contents  of  the  column  were  placed  into  sealed 
plastic  tubes.  The  aqueous  solution  in  the  tubes  was  poured  into  a  sample  vile  and 
weighed.  1  mL  of  nitric  acid  was  added  to  9  mL  of  this  solution.  EPA  method  3015 
(USEPA  1997)  was  then  followed  to  find  the  concentration  of  Ag  in  the  solution. 


20 


Ten  to  15  mL  of  nitric  acid  was  added  to  the  sand  from  the  tubes,  and  the  samples 
where  centrifuged  for  10  minutes  at  2000  rpm.  The  supernatant  was  poured  from  the 
sample  and  weighed.  lmL  of  nitric  acid  was  then  added  to  9  mL  of  supernatant.  These 
samples  also  were  analyzed  using  EPA  method  3015  (USEPA  1997)  on  the  ICP. 

Finally,  to  find  the  amount  of  AgNPs  attached  to  the  sand  approximately  1  g  of 
sand  was  taken  from  each  sample  after  the  supernatant  was  removed  and  it  was  placed  on 
a  watch  glass  and  placed  in  an  oven  at  105  °C  for  12  hours.  Then  0.5  g  of  the  dried  sand 
was  put  into  a  sample  container  with  4.5  mL  of  HNO3,  and  0.5  mL  hydrochloric  acid. 
These  samples  were  analyzed  with  EPA  method  3051  (USEPA  1997). 

Mass  Balance  Experiment 

An  additional  experiment  was  run  using  coarse  sand  and  a  4  mL  min'1  flow  rate  to 
evaluate  the  overall  silver  mass  balance.  The  column  experiment  was  set  up  and  run  as 
described  above;  however  in  the  previous  experiments  once  the  effluent  of  the  column 
passed  thru  the  UV-vis  for  measurement  of  the  AgNPs  concentration,  it  was  weighed  and 
then  discarded.  For  this  experiment,  all  90  effluent  samples  were  retained  for  further 
analysis.  The  odd  sample  numbers  were  prepared  for  inductively  coupled  argon  plasma- 
optical  emission  spectroscopy  (ICP-OES)  analysis.  Each  of  these  samples  contained 
approximately  2  mL  of  solution.  From  this,  1  mL  was  taken  and  9  mL  of  DI  water  was 
added.  Then  lmL  of  HNO3  was  added  to  9  mL  of  the  diluted  solution.  The  samples 
were  then  microwave  digested  using  EPA  method  3015  (USEPA  1997)  and  analyzed  on 
the  ICP-OES.  This  method  would  account  for  both  dissolved  and  particulate  silver,  in 
contrast  to  the  UV-vis  measurement,  which  would  only  account  for  particulate  silver. 
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Influent  samples,  nominally  at  a  concentration  Co  =  15  mg  L  ,  were  also 
analyzed  using  the  ICP.  After  the  experiment  was  run,  9  mL  of  the  remaining  influent 
solution  was  taken  and  lmL  of  HNO3  was  added.  After  microwave  digestion,  1  mL  was 
removed  and  added  to  9  mL  of  RO  water  to  achieve  a  dilution  factor  of  10.  The  sample 
was  then  run  on  the  ICP-OES,  to  quantify  Co. 

The  sand  in  this  mass  balance  experiment  was  analyzed  as  in  the  previous 
experiments  with  one  addition.  For  this  experiment,  the  sand  was  centrifuged  an 
additional  time  with  15  mL  of  nitric  acid  for  10  minutes  at  2000rpm.  Again,  the 
supernatant  was  poured  from  the  sample  and  weighed.  lmL  of  nitric  acid  was  then 
added  to  9  mL  of  supernatant.  This  was  accomplished  to  ensure  all  the  captured  silver  in 
the  sand  would  be  released  to  be  quantified. 

Modeling 

To  help  us  better  describe  and  understand  AgNP  transport  in  porous  media,  a 
model  was  developed.  The  model  incorporates  the  processes  of  advection,  dispersion, 
and  capture  to  describe  transport  in  a  porous  medium.  AgNP  capture  was  modeled  as 
both  reversible  (2nd  term  on  the  left-hand  side  of  equation  1)  and  irreversible  (third  tenn 
on  the  right-hand  side  of  Equation  1) 


adN  dN s  nd2N  ndN 

9 —  +  p — -  =  DO — —~v0 - AON 

8t  dt  dx  dx 


Eq.  1 


Reversible  capture  was  modeled  as  either  a  fast  equilibrium  process  (Equation  2) 
or  a  rate-limited  process  (Equation  3). 


a 


Eq.  2 
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dN  kf 

—^  =  a(^N-Ns)  Eq.3 

dt  a 

where  N  and  Ns  are  the  concentrations  of  suspended  and  captured  nanoparticles, 
respectively,  p  is  the  bulk  density  of  the  porous  medium,  0  is  the  porosity  of  the  porous 
medium,  D  is  a  dispersion  coefficient,  v  is  the  average  pore  velocity  of  the  water,  t  is 
time,  x  is  length,  X  is  a  first-order  irreversible  capture  rate  constant,  and  kf  and  a  are  first- 
order  capture/release  rate  constants,  respectively.  For  a  clean  one-dimensional  column 
with  a  finite-pulse  third  type  boundary  condition  at  the  inlet,  the  initial  and  boundary 
conditions  are: 


N(x,t  =  0)  =  Ns(x,t  =  0)  =  0 

(vAT-D^)!  =vN  0<t<t 

OX 

=  0  t>tp 


Eq.  4a 

Eq.  4b 


N(x  — »  oo,  f)  =  0 


Eq  .  4c 


where  nanoparticles  at  a  suspended  concentration  of  No  are  injected  into  the  column  inlet 
(x  =  0)  for  a  time  period  tp. 

A  nondimensional  version  of  Equations  1  and  2  is 


R 


dN  _  1  d2N  dN 
dT  Pe  dx2  dx 


- Da 1  N 

irrev 


Eq.  5 


where 

T  =  —  is  a  dimensionless  time  scale,  also  referred  to  as  pore  volume  (PV) 


x 


x 

L 


is  a  dimensionless  length  scale 
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L  is  the  length  of  the  column 
—  N 

N  = - is  a  nondimensional  suspended  nanoparticle  concentration 

N0 

pn 

R  =  1  +  — —  is  a  retardation  factor 
9 


Pe  =  —  is  a  Peclet  number,  the  ratio  of  a  dispersion  time  scale  to  an  advection 
time  scale 

/LL 

Dctlrev  =  —  is  a  Damkohler  number,  the  ratio  of  an  advection  time  scale  to  an 
v 

irreversible  capture  time  scale 

Nondimensional  versions  of  Equations  1  and  3  are 


dN  dN s  1  d2N  dN  n  7  - 

- 1 - — - ~ - Da  :  N 

ST  GT  Pe  ox2  dx 


Eq.  6 


dN 
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Eq.  7 


where 


aL  . 


Darev  =  — is  also  a  Damkohler  number,  in  this  case  the  ratio  of  an  advection 
v 


time  scale  to  a  reversible  capture/release  time  scale 


-  N 
N..  = 


S  k 


N  n 

a  0 


is  a  nondimensional  attached  nanoparticle  concentration 


The  nondimensionalized  initial  and  boundary  conditions  (corresponding  to 
Equation  4)  are: 
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N(x,T  =  0)  =  Ns(x,T  =  0)  =  0 


Eq.  8a 


(A,-^^)|,.o=l  0<T<T 

Pe  ox 

=  0  T>Tp 
N(x  — >  go,  T)  =  0 


Eq.  8b 

Eq  .  8c 


where  Tp  =  — —  is  the  dimensionless  input  pulse  time 
L 

For  all  the  UV-Vis  results,  we  use  an  analytical  solution  to  Equations  5  and  8  to  plot  N 
vs  T.  Pe,  R,  P>a\rrev ,  and  Tp  are  used  as  fitting  parameters. 

For  the  ICP  results,  we  use  an  analytical  solution  to  Equations  6,  7,  and  8 
(Valocchi  and  Werth,  2004)  to  plot  N  vs  T.  Pe,  Da'rev ,  Da\rrev ,  and  Tp  are  used  as  fitting 
parameters. 


Results  and  Discussion 

Characterization  of  AgNPs 

The  HRTEM  image  and  analysis  (Figures  la  and  lb)  show  results  similar  to  those 
presented  by  Flory  et  al  (2013).  At  the  neutral  pH  used  in  this  study,  the  diameters  of  all 
particles  ranged  from  4  to  38  nm  with  95%  of  the  particles  under  25  nm,  4%  between  25- 
30  nm  and  less  than  1%  over  3 1  nm.  The  average  diameter  was  16.9  nm. 
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Characterization  of  Sand 

XRD  showed  that  all  sands  used  were  similar  to  quartz.  No  metal  oxides  were 
found  in  the  crystalline  form  attached  to  any  of  the  sand  size  fractions.  Amorphous  phase 
oxides  could  be  present. 

X-ray  photoelectron  spectroscopy  (XPS)  was  performed  on  the  AgNPs  to  confirm 
the  chemical  composition  of  the  nanoparticles ’  surface.  It  is  found  that  the  top  three 
nanometers  of  the  surface  were  composed  of  75%  carbon,  13%  oxygen,  10%  nitrogen 
and  1%  metallic  Ag.  The  presence  of  carbon  is  a  result  of  the  coating  applied  to  the 
AgNP  by  the  manufacturer. 

Transport  of  AgNPs 

The  dimensionless  breakthrough  curves  for  the  chloride  tracer  (Figure  2a)  through 
the  coarse  sand  column  at  1  and  4  mL  min'1  are  shown.  Shown  in  Figure  2b  and  2c  are 
the  AgNP’s  BTCs  at  the  1  and  4  mL  min'1  flow  rate  through  the  coarse  sand.  Figure  3a 
plots  the  BTCs  for  the  chloride  tracer  through  the  medium  sand  column  at  1  and  4  mL 
min'1.  Figures  3b  and  3c  show  the  BTCs  for  AgNPs  through  the  medium  sand  column  at 
1  and  4  mL  min'1.  Table  1  shows  the  zeroth  and  first  moment  for  the  BTCs.  No 
breakthrough  of  AgNPs  was  observed  for  the  fine  sand  at  either  flow  rate.  As  a  result,  no 
zeroth  or  first  moment  is  shown  for  the  fine  sand. 

In  Table  2  the  ratio  of  the  zeroth  and  first  moments  of  the  AgNP  BTCs  to  the 
moment  of  the  corresponding  tracer  BTC  is  presented.  The  zeroth  moment  ratio  shows 
the  extent  of  capture  of  the  AgNPs  compared  to  the  tracer,  which  is  assumed  to  not  be 
captured  at  all.  The  first  moment  ratio  shows  the  extent  of  retardation  compared  to  the 


26 


tracer,  which  is  also  assumed  to  not  be  retarded  at  all.  Comparing  the  ratio  of  the 
moments  at  the  two  flow  rates,  it  can  be  seen  that  at  least  for  the  fourfold  increase  in  flow 
studied  here,  the  flow  rate  had  little  effect  on  the  transport  of  AgNPs.  In  addition,  the 
retardation  did  not  change  significantly  with  the  different  flow  rates  in  the  different  types 
of  sand.  The  retardation  observed  here  differed  from  what  Flory  et  al.  (2013)  observed. 
Flory  et  al.  (2013)  found  no  AgNP  retardation  in  their  glass  bead  media  whereas  some 
retardation  of  the  AgNPs  was  observed  in  the  current  study  with  different  types  of  sand. 

When  comparing  the  effects  of  media  size,  the  medium  sand  (0.21-0.297  mm)  had 
less  retention  of  AgNPs  than  the  coarse  sand  (0.4-2  mm)  in  these  experiments  as  shown 
by  the  zeroth  moment  ratios  in  Table  2  as  well  as  the  data  in  Figure  4.  This  is  not  in 
accordance  with  conventional  filtration  theory  which  states  that  with  smaller  grain  size 
the  fraction  of  mass  retained  in  the  media  should  increase.  This  also  disagrees  with  Li  et 
al.  (2008)  who  found  that  retention  was  increased  as  the  particle  size  of  the  media 
decreased.  We  speculate  that  our  observation  may  be  attributed  to  the  difference  in  the 
size  range  of  the  medium  vs  coarse  sand,  as  opposed  to  the  median  grain  size  difference. 
The  coarse  sand  had  a  size  range  from  0.4-2  mm,  and  the  medium  sand  had  a  much 
narrower  range  of  0.21-0.291  mm.  Zhang  et  al.  (2012)  showed  similar  results  in  their 
transport  studies  with  Fullerene  nanoparticles  in  Lula  soil.  They  found  that  significantly 
more  mass  was  retained  in  Lula  soil  versus  Ottawa  sand  (mainly  pure  quartz).  They 
attribute  this  to  the  irregular  and  rougher  shape  and  greater  heterogeneity  of  Lula  sand 
compared  with  Ottawa  sand. 
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Ren  and  Smith  (2013)  concluded  the  size  of  the  media  had  an  effect  on  the 
transport  of  AgNPs.  They  found  increased  retention  with  smaller  media  diameters.  In 
their  study  their  coarse  and  medium  sand  median  grain  diameters  (0.55  and  0.34  mm) 
were  both  within  the  range  of  our  coarse  sand  (0.2-4  mm).  The  fine  sand  used  in  our 
experiment  had  a  diameter  of  <  0.074  mm.  The  result  obtained  from  our  fine  sand 
experiments  agree  with  filtration  theory  as  well  as  with  Li  et  al.  (2008)  and  Ren  and 
Smith  (2013). 

As  in  the  other  experiments,  a  BTC  was  developed  for  the  single  mass  balance 
experiment.  We  found  that  17%  of  the  initial  mass  passed  through  the  column.  We  also 
found  7 1 .5%  of  the  silver  was  captured  by  the  sand.  A  total  of  88.75%  of  the  mass  was 
accounted  for  in  this  experiment. 

UV-vis  measurements  were  taken  along  with  the  ICP-OES  measurements  (Figure 
5).  It  is  clear  that  more  mass  is  measured  in  the  effluent  with  the  ICP-OES.  The 
difference  in  mass  can  be  attributed  to  the  UV-vis  only  measuring  AgNPs  rather  than  the 
total  silver  as  the  ICP-OES  does.  It  appears  that  much  of  the  mass  being  detected  is  in 
the  dissolved  form.  Also  the  ICP-OES  data  were  able  to  show  that  tailing  occurred  as  the 
experiment  continued.  This  tailing  effect  can  be  attributed  to  the  captured  silver  being 
dissolved  and  then  measured  by  the  ICP.  This  dissolved  silver  tail  is  not  detected  by  the 
UV-vis,  which  only  detects  particulates. 

It  was  noted  that  over  the  several-month  course  of  these  experiments,  AgNP 
breakthrough  curves  could  not  be  replicated,  even  though  experimental  conditions  were 
ostensibly  the  same.  Experiments  were  run  five  months  after  the  initial  experiments  in  an 
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attempt  to  reproduce  results  previously  obtained.  It  can  be  seen  for  the  coarse  sand  at  4 
mL  min'1  flow  rate  column  experiments  run  five  months  apart,  there  was  much  less 
breakthrough  of  AgNPs  measured  by  the  UV-vis  after  five  months  (Figure  6).  The 
samples  taken  at  the  later  experiment  were  also  analyzed  by  ICP-OES  with 
approximately  89%  of  the  mass  accounted  for.  The  influent  concentration  measured  by 
UV-vis  for  this  experiment  was  similar  to  all  other  experiments  run  (-15  mgL'1); 
however  it  is  clear  that  the  effluent  concentrations  are  not  similar.  One  hypothesis  is  that 
the  PVP  coating  for  the  AgNPs  has  deteriorated  after  five  months  and  made  the  particles 
easier  to  dissolve  into  solution.  Another  hypothesis  could  be  that  the  PVP  deteriorated, 
causing  the  NPs  to  agglomerate  and  be  captured  by  the  media.  Then,  the  captured 
particles  dissolved,  resulting  in  the  long  tail  observed  in  the  ICP.  These  hypotheses  could 
not  be  tested  at  this  time. 

Model  Simulations  of  AgNP  Transport 

As  can  be  seen  in  Figures  2-4,  the  equilibrium  model  does  a  good  job  in 
simulating  the  experimental  data.  We  also  see,  by  comparing  Tables  2  and  3,  that,  as 
should  be  expected,  the  Damkohler  number,  which  indicates  irreversible  capture,  is  larger 
when  the  0th  moment  ratio  is  smaller,  and  the  retardation  factor  is  larger  when  the  first 
moment  ratio  is  larger. 

From  Figure  7,  we  see  that  to  simulate  the  transport  of  total  silver,  it  is  necessary 
to  apply  a  model  that  assumes  both  irreversible  capture  and  rate-limited  reversible 
capture.  Figure  7  plots  the  BTC  for  the  coarse  4  mL  min'1  flow  rate  and  coarse  sand  with 
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data  obtained  from  the  ICP-OES  analysis.  This  graph  shows  that  there  is  tailing  of  the 
silver  in  the  experiments.  This  tailing  was  not  observed  in  the  other  experiments,  which 
used  the  UV-vis  to  measure  silver,  because  the  limits  of  detection  for  the  UV-vis  were 
higher  than  those  of  the  ICP  and  also  the  UV-vis  only  detects  particles  whereas  the  ICP 
measures  total  silver. 

Significance 

The  effect  of  flow  rate  and  media  size  on  the  transport  of  AgNPs  under  saturated 
conditions  was  examined  in  this  study.  It  was  found  through  experimentation  and 
mathematical  modeling  that  the  flow  rate  plays  a  minor  role  in  transport  compared  to  the 
size  of  the  media.  More  specifically,  the  size  unifonnity  of  the  media  appears  to  play  the 
larger  role  in  affecting  transport. 

We  found  that  while  a  medium  can  have  larger  average  diameters  of  particles,  it 
can  retard  and  capture  more  AgNPs  than  a  finer  sand  that  is  better  sorted.  This  is  relevant 
in  the  natural  environment  because  of  the  differences  in  media  that  may  be  present  at 
different  locations.  We  also  confirmed  previous  studies  that  AgNPs  are  filtered  out  by 
fine  media  (<0.074  mm). 
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III.  Conclusion 


Summary  of  Findings 

The  effect  of  flow  rate  and  media  size  on  the  transport  of  AgNPs  under  saturated 
conditions  was  examined  in  this  study.  It  was  found  through  experimentation  and 
mathematical  modeling  that  the  flow  rate  plays  a  minor  role  in  transport  compared  to  the 
size  of  the  media.  More  specifically,  the  size  uniformity  of  the  media  appears  to  play  the 
larger  role  in  affecting  transport. 

We  found  that  while  a  medium  can  have  larger  average  diameters  of  particles,  it 
can  retard  and  capture  more  AgNPs  than  a  finer  sand  that  is  better  sorted.  This  is  relevant 
in  the  natural  environment  because  of  the  differences  in  media  that  may  be  present  at 
different  locations.  We  also  confirmed  previous  studies  that  AgNPs  are  filtered  out  by 
fine  media  (<0.074  mm). 

Recommendations  for  Future  Research 

Due  to  the  numerous  physical,  chemical,  and  even  biological  factors  that  may 
affect  transport  of  ENPs  in  the  environment,  along  with  the  rapidly  expanding  uses  of 
these  materials,  it  is  important  that  we  gain  a  better  understanding  of  how  these  factors 
affect  ENP  fate  and  transport.  While  not  conclusive,  observations  in  this  study  point  to 
media  heterogeneity  as  an  important  factor  which  affects  transport  of  AgNPs.  To  test  this 
hypothesis,  an  experiment  could  be  conducted  whereby  the  coarse  sand  used  in  this  study 
is  sorted  into  a  more  homogenous  sand,  and  AgNP  transport  in  this  more  homogeneous 
media  compared  to  transport  in  the  unsorted  heterogeneous  media.  Key  to  understanding 
would  be  quantifying  the  relative  heterogeneities  of  the  different  media,  perhaps  through 
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use  of  a  uniformity  coefficient.  Various  tests  could  then  be  run,  comparing  transport  of 
ENPs  through  media  with  the  same  median  diameters,  but  with  different  coefficients  of 
uniformity. 

It  was  noted  that  the  AgNP  breakthrough  behavior  toward  the  end  of  the  study 
was  significantly  different  than  the  behavior  observed  months  earlier,  at  the  study’s  start, 
even  though  conditions  (media  size,  flow  rate)  were  identical.  It  was  hypothesized  that 
over  the  months  the  AgNP’s  PVP  coating  had  changed  in  some  manner  as  to  affect  their 
transport.  Experiments  to  detennine  if  and  how  the  AgNP  characteristics  change  over 
time,  and  how  these  changes  affect  breakthrough  behavior,  would  be  very  useful  in 
helping  us  to  understand  how  these  ENPs  are  transported  in  the  environment. 
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Table  1 :  Breakthrough  curve  moments  for  AgNP  at  different  flow  rates  and  sand  sizes 


0th  moment 

1st  moment 

Coarse  4  mL  min'1  tracer 

4.64 

4.07 

Coarse  4  mL  min'1 

2.14 

4.83 

Coarse  1  mL  min"1  tracer 

4.11 

4.25 

Coarse  1  mL  min"1 

1.69 

5.11 

Med  4  mL  min' ‘tracer 

4.75 

4.84 

Med  4  mL  min'1 

2.79 

5.25 

Med  1  mL  min'1  tracer 

5.41 

4.75 

Med  1  mL  min'1 

3.28 

4.84 

Table  2:  Ratio  of  AgNP  to  tracer  breakthrough  curve  moments  at  different  flow  rates  and 
sand  sizes 

0th  moment  ratio 

1st  moment  ratio 

Coarse  4  mL  min"1 

0.46 

1.19 

Coarse  1  mL  min"1 

0.41 

1.20 

Med  4  mL  min'1 

0.58 

1.08 

Med  1  mL  min'1 

0.61 

1.02 

Table  3:  Model  parameter  values  at  different  flow  rates  and  sand  sizes 

R 

Pe 

rrev 

Tp 

Coarse  1  mL  min"1 

3.30 

5.26 

0.97 

4.33 

Coarse  4  mL  min’1 

2.64 

6.67 

0.77 

4.63 

Medium  1  mL  min'1 

1.15 

0.46 

0.40 

6.9 

Medium  4  mL  min’1 

2.43 

1.92 

0.71 

6.25 
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Figure  1 :  a)  HRTEM  images  of  AgNPs  and  b)  size  histograms  prepared  based  on 
analysis  of  the  HRTEM  images  of  >  500  AgNPs. 
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a)  1 


C/Co 


♦  1  mL/min 

Model  1  mL/min 


♦  4  mL/min 

Model  4  mL/min 


Figure  2:  Coarse  sand  BTC  for  a)  tracer  data  and  model  simulation  for  1  and  4  mL  min'1 
flow  rates  b)  AgNP  data  and  model  simulation  for  1  mL  min'1  flow  rate  c)  AgNP  data 
and  model  simulation  for  4  mL  min'1  flow  rate 
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Figure  3:  Medium  sand  BTC  for  a)  tracer  data  and  model  simulation  for  1  and  4  mL  min'1 
flow  rates  b)  AgNP  data  and  model  simulation  for  1  mL  min'1  flow  rate  c)  AgNP  data 
and  model  simulation  for  4  mL  min'1  flow  rate 
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Figure  4:  BTC  for  AgNPs  through  coarse  and  medium  sand  columns  and  1  mL  min'1 


flow  rates  with  the  corresponding  model  simualtions 
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Figure  5:  AgNP  BTCs  for  coarse  sand  column  and  4  mL  min'1  flow  rates  plotted 
comparing  ICP  and  UV-vis  data 
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Figure  6:  AgNP  BTCs  for  the  coarse  sand  column  and  a  4  mL  min'1  flow  rate  taken  at 
time  zero  and  after  five  months  measured  with  UV-Vis 


Figure  7:  AgNP  BTC  data  measured  on  the  ICP-OES  and  model  simulation  for  the 
coarse  sand  column  and  a  4  mL  min’1  flow  rate 
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